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Adiabaticity of quantum evolution is important in many settings; one example is the adiabatic quantum com- 
putation (AQC). Nevertheless, to date, there is no effective method available for testing the adiabaticity of the 
evolution for the case where the eigenenergies of the driven Hamiltonian are not known. We propose a simple 
method for checking adiabaticity of a quantum process for an arbitrary quantum system. We further propose an 
operational method for finding more efficient protocols that approximate adiabaticity, and suggest a "uniformly 
adiabatic" quench scheme based on the Kibble-Zurek mechanism for the case where the initial and the final 
Hamiltonians are given. This method should help in implementing AQC and other tasks where preserving the 
system in the ground state of a time-dependent Hamiltonian is desired. 
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In a quantum quench process, when the Hamiltonian of a quantum system is driven from Hq to Hi, interstate excitations of 
r - — , , the system usually occur, owing to the non-commutativity of the Hamiltonians at different moments. However, when the quench 
^G ' process is slow enough, the interstate excitations will be suppressed. According to the quantum adiabatic theorem yj], when the 

ipi condition for quantum adiabatic approximation ( Aground ft) I d \ ■. \ ^excited (t) \ <C A(i) 2 is satisfied, the system will remain in 
r* , the ground state - its evolution will be adiabatic - except for some special situations yj. Here H(t) is the changing Hamiltonian, 
and A(i) is the minimal energy gap between the ground state |$ gr ound(i)) an( ^ tne fi rst excited state | < P xcited(i)) of H (t). 

In order to ensure that a quantum system evolves adiabatically, one usually needs to find the energy spectrum (or at least the 
smallest energy gap A) of the driven Hamiltonian. One can then use the quantum adiabatic theorem to choose a proper time 
scale, so that the conditions for quantum adiabatic approximation are satisfied and the evolution remains adiabatic. Nevertheless, 
C*") ■ in practice, neither eigenenergies nor eigenstates of a complex quantum many-body system are easy to obtain. This is often the 
case in implementing the adiabatic quantum computation (AQC) la] as well as quantum annealing 0,1a]- Hence, one does not 
have the ingredients to use quantum adiabatic theorem. One cannot count on the direct comparison between the final state and 
the instantaneous ground state of the final Hamiltonian, either. Thus, it would be useful to find a reliable method for evaluating 
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I. INTRODUCTION 



J^-: , the adiabaticity of an evolution under an arbitrary Hamiltonian, especially when one has no idea about the eigenstates and/or 
eigenenergies of the system (except at the initial moment). 

The quench echo method we propose here is one solution to the above problem. It will allow one to evaluate unambiguously 
the adiabaticity of a process. What is more, it can help one find the efficient "uniformly adiabatic" quench path in the parameter 
space of the Hamiltonian. Such ideas may have applications in the implementation of AQC jg]. This paper is organized as 
follows: In section II, we introduce the quench echo method and briefly explain its underlying physics. In section III, we use a 
simple model to demonstrate main ideas of the general theory. In Section IV, we propose a uniformly adiabatic scheme that is 
based on the application of the Kibble-Zurek mechanism (KZM) to quantum phase transitions. In Section V we give discussions 
lJ ■ and conclusions. 

II. QUENCH ECHO 

Consider a system described by the Hamiltonian H(g(t)), where g(t) is a time-dependent parameter. The system is initially 
prepared in the ground state of H(g(t = 0)). The system evolves under the influence of the driven Hamiltonian, which changes 
from H(g(t = 0)) to H(g(t = T)). Our aim is to test the adiabaticity of this evolution, but we know nothing about the 
eigenenergies and eigenstates of the time-dependent Hamiltonian except at t — 0. Hence we cannot count on the comparison 
between the evolving state and the instantaneous ground state. Neither can we use the adiabatic theorem. Nevertheless, we can 
apply a backward "echo" quench following the initial quench (from t = to t = T). That is, one extends the evolution from 
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t = Ttot = 2T$$\: 



, H(g(t)),(0<t<T) 

H ^>>-< H(g(2T-t)),(T<t<2T) ' (1) 



The final Hamiltonian is identical to the initial Hamiltonian H(g(t = 2T)) = H(g(t = 0)). Hence, we can use the fidelity of the 
initial state, e.g., the ground state of H(g(t = 0)) and the final evolving state as a criterion for the adiabaticity of the evolution. 
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where f is the time-ordered operator and \GS) is the ground state of H{g(t = 0)). When the fidelity F is greater than a threshold 
value close to unity, (e.g., 0.999, the error tolerance is 0.001), the whole process (0 < t < 2T) is adiabatic. This implies that the 
forward quench process (0 < t < T) is adiabatic. The underlying mechanism for this "quench echo" method is straightforward: 
Except for the eigenstate of the Hamiltonian at the initial moment H(g(t = 0)), we do not have any information about its 
eigenstates at other moments. Hence, we can only quench the Hamiltonian back, so that it goes back to its initial H(g(t = 0)), 
and we can measure the final state (and compare it with the initial state). The quench echo protocol (Eq. dTJ) ensures that the 
excitation probabilities in the forward quench process and those in the backward quench process are similar (but not identical; see 
Refs. IUll[l2ll ). As a result, when the forward process is adiabatic (no excitations), so is the backward quench process. Otherwise 
both the forward and the backward processes are nonadiabatic, and the phase accumulated between transitions (known as the 
Stiickelberg phase) may result in constructive or destructive interference ||5l UM- Usually the excitations in the forward and 
the backward processes cannot cancel each other out (but see Ref. ||5|]). Hence, through this quench echo method, without 
knowing about the eigenenergies and eigenstates, one is able to evaluate the adiabaticity of an arbitrary evolution in most cases. 
Nevertheless, in some special cases (e. g., impulse evolution, which is so fast that the state of the system is frozen) the final 
fidelity is equal to unity, but the process may not be adiabatic. A solution to this problem is to let the system evolve freely for 
some time before the backward quench. We will discuss this in detail in the next section. By utilizing quench echo one can even 
find a uniformly adiabatic quench protocol for a given Hamiltonian by repeating the above process with different quench time 
scales. 

III. A CASE STUDY: ISING CHAIN IN A TRANSVERSE MAGNETIC FIELD 

We now use a simple model to demonstrate our central ideas: Consider the quench dynamics of an Ising model in an transverse 
magnetic field 114JI . The time-dependent Hamiltonian is 

N 

H(t) = -jJ2\v?v! + i+9(t)^}, (3) 

where J indicates the energy scale; erf, a — x, y, z is the Pauli matrix on the i th lattice site; and g(i) is the reduced strength of 
the magnetic field, which varies with time. It is known that for this model there is a finite energy gap A = 2 Jjj at g — cos(w/N) 
when the size of the system N is finite. For simplicity, we consider a linear quench protocol 

{.go , (0 < t < (.go - 9t)tq), 
Zgr-ga H , ((30 - 9t)tq < t < 2(g Q - g T )T Q ), 
T Q 

where tq is the time scale of the quench. The larger the tq, the slower the quench. In the forward quench the strength of the 
magnetic field is ramped from g = g to g = g? continuously, and in the quench echo, it is ramped back from g = g^ to g = go- 
where gr is the turnaround point. Initially the system is prepared in the ground state of H(g = go). When one quenches the 
system at different rates (by choosing different tq), the fidelity $2& will be different. 

The Hamiltonian of the Ising model (0 can be decoupled into N independent fermionic modes 11411 . 

H(t) = 5>*(fl(t)) [| + (i)>* (+(t)| fe - \-(t)) k (-(i)U , (5) 

fc 

where \ + {t)) k and | — (t)) k are the two instantaneous eigenstates of the k mode. Their corresponding eigenenergies are 
±A k (g(t)), and A{g(t)) = J^/g{t) 2 - 2g(t) cosfc + 1. Here k = (2s + 1)tt/A, s = 0,1,2,- •• ,N/2 - 1 is the wave 
vector, and the number of spins N is even. 
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FIG. 1: Three regimes of the quench dynamics: (a) nearly adiabatic regime, tq = 150; (b) the intermediate regime tq — 35; and (c) nearly 
impulse regime tq — 0.004. The horizontal axis represents the parameter <?(£), which varies between and 10, and the vertical axis represents 
the probability of the system being in the instantaneous ground state during the evolution. There is a quantum phase transition at g c = 1. The 
red solid line represents the forward quench (from <?n = 10 to gr = 0) and the green dashed line represents the quench echo (from gr = 
to go — 10). In both the nearly adiabatic and the nearly impulse regimes, the fidelity at final moment is close to unity. In the first row (a)-(c), 
there is no time delay at the turnaround point. In the second row (a')-(c'), the delay time at gr = is At — 0.1, 0.2, 0.3, 0.4. In the third row 
(a")-(c"), the delay time at gr = is At = 10, 20, 30, 40. The number of the spins in the Ising chain is iV = 50. 

We write the Schrodinger equation ih-S^ |$(i)) = H{t) |$(i)) in the instantaneous eigenbases \+(t)) k and | — (t)) k of H(t), 
where |$(i)) = Yik a k(t) !+(*))& + Pk(t) \— (£)}&• F° r simplicity we choose ft = 1 hereafter. For both the forward quench 
(0 < t < (.go — 9t)tq) and the backward ((go — cjt)tq < t < 2(go — 9t)tq) process , the Schrodinger equation can be written 

as 



dt 



a k (t) 
Pk(t) 



2A fc (ff(t)), 



-iJ 3 sink dg(i) 



2kt(g(t)) dt 



i J 3 sink dg(t) 
2A 2 Mt)) dt 



,-2A fe Q?(£)) 



afc(t) 

Pk{t) 



(6) 



where the initial condition for Eq. (O is ak(t = 0) = 0, /3k(t = 0) = 1. The modulus square of the overlap between the final 
state of Eq. (O and the instantaneous ground state at g = go gives the fidelity (O 



F = P G 5(2(,9o - g T )T Q ) = \\ W(g Q - 9 t)tq)\' 



(7) 



fc>0 



In the following, we will focus on the solution of the Eq. ©. We will consider both the numerical and the analytical results. 



A. Kibble-Zurek mechanism and three regimes 

Before the quantitative study of the fidelity and its relation with the time scales of the quench, we describe the Kibble- 
Zurek mechanism (KZM) lfl5l[l6ll of second-order phase transitions, which provides a quantitative understanding of the quench 
process. The KZM describes e.g., the relation between the density of topological defects, which are generated during quenching 
across a phase transition, and the time scale of the quench (see Ref. IU8I1 for a good review). The KZM was first introduced in 
the classical phase transitions lfl5l[l6ll . and later generalized to quantum phase transitions 12111 . In our study, however, we will 
not focus on the density of topological defects, but on the adiabaticity of the evolution of the system. 

A quantum phase transition is characterized by a vanishing excitation gap A(g(t)) ps Ao\g(t) — g c \ zv and a divergent 
correlation length £ sa £o/ls(£) — 9c\ v , where z and v are the critical exponents, and Ao and £o are constants IU4I1 . We define a 
dimensionless distance from the critical point g c by 

e(t) = iVh^L. (8) 



A generic e(t) can be linearized near the critical point e(t) = as 01711 : 

e(t)«- — . (9) 

There are two interlinked time scales during a quench: the system reaction time given by the inverse of the gap r(e(i)) = 
1/A \g{t) — g c \ zv and the time scale of transition given by \g(t) — g c \ z,y /jjj\g(t) — g c \ zu ■ Away from the critical point the 
reaction time is small in comparison with the time scale of transition and the evolution is adiabatic. Near the critical point, 
however, the opposite situation occurs and the evolution is approximately impulse (the state of the system is frozen out). The 
boundary t between the two regions is determined by the relation r(e(t)) = e/e| 4 ~, or 

i _ M-9.\?L m (10) 



\9{t)-9c\ zv i\g(i)-g c \ z »' 

That is, i ~ (-f-) zu , which gives i ~ Tq + " 1 ' ITmulll . For the Ising model, we have z = v = 1 resulting in i ~ Tq lulTl . 
According to KZM when t <G (— i, i), the system will not evolve (the wavefunction will be frozen). Outside this time interval 
the system will evolve approximately adiabatically. 

For an infinitely large system, the energy gap is vanishingly small at the critical point, and one can always find a t. According 
to the KZM this implies that, no matter how slow one quenches the Hamiltonian in an infinite system, the evolution across the 
critical point can never be adiabatic. For a finite-size system, however, there is a finite energy gap even at the critical point. When 
one quenches the system sufficiently slowly (large tq), t approaches very near the critical point where - for a finite systems - 
scalings no longer hold. As a consequence, the KZM does not lead to simple scaling, as r(e(i)) = e/e|g leads to a more difficult 
equation which has to be solved to obtain i. 12211 . Indeed - in accord with the adiabatic theorem - the KZM predicts that when tq 
is larger than the inverse of the gap the transition will remain adiabatic throughout. Thus, a finite energy gap allows an adiabatic 
evolution across the critical point when the Hamiltonian is driven sufficiently slowly. This is the adiabatic quench regime. By 
contrast, when one quenches the Hamiltonian very fast (small tq), there is a big t and there is an approximately impulse regime 
for t E (— t, i) when the quench is essentially instantaneous. In this time interval the system will approximately cease to evolve 
- its wavefunction will be frozen. This is the so-called impulse regime Il2ll 12311 . When one chooses a time scale of quench tq 
between the above two limiting cases, the system will evolve adiabatically when either t < — t or t > t, and will be frozen when 
t € (— t, i). We call this regime the intermediate regime. We can summarize the quench behavior as follows: For a finite-size 
system, when tq is large enough, the evolution will be adiabatic; When tq is extremely small, the state of the system will be 
frozen; When tq is in between these two limiting cases, the process is in the intermediate regime. 

B. Numerical and analytical results 

Having obtained the qualitative understanding of the quench dynamics from the above KZM arguments, in the following we 
will study the Ising model quantitatively, and compare the results with the estimates obtained above. We consider a spin chain 
with a finite size N = 50, and start evolving it at go = 10 and let it turn around at gT = 0. There is a finite energy gap for this 
system at the quantum critical point g c = 1. We choose three different quench time scales tq = 150, tq — 35, and tq = 0.004 
which correspond to the adiabatic, intermediate, and impulse regimes. The system evolves under the time-dependent Schrodinger 
equation. We plot the probability Pos(g) in the instantaneous ground state as a function of the controlling parameter g during 
the quench process in Fig. la-lc. 
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FIG. 2: The fidelity as a function of the time scale of quench tq. Panels (a)-(d) represent different free evolution time at turnaround point 
gr — before the quench echo. Here N = 50, and start point at go = 10, and the times of free evolutions are chosen to be At = 0, At = 0.1, 
At = 0.3, At = 0.7. It can be seen that in the adiabatic regime the fidelity is always equal to unity, but in the impulse regime it is less than 
unity and varies with the time of free evolution At. 



From Fig. la it can be seen that when the time scale of the quench is relatively large, the system evolves almost adiabatically 
in the whole range of the parameter gx = < g < go = 10, Pcs(g) is always close to unity (except a tiny decay and partial 
revival at the critical points. So is the fidelity of the quench echo (see Fig. la). 

When the time scale of the quench is reduced to tq = 35 (see Fig. lb), the quench dynamics enters the intermediate regime. 
It can be seen that away from the critical point, the evolution is adiabatic. But near the critical point, the probability in the 
instantaneous ground state Pcs(g) decays sharply and oscillates rapidly. This is due to the interstate transitions at the anti-cross 
point. Soon after passing through the quantum critical point the adiabatic evolution resumes. 

When the time scale of the quench is further reduced to tq = 0.004 (almost instantaneous quench), the wave function of the 
system is nearly frozen. Hence, the probability of being in the instantaneous ground state is simply equal to the overlap of the 
initial state and the instantaneous ground state. In the backward quench, the same situation arises. Because in both the forward 
and the backward quench, the wave functions of the system are frozen, and hence are almost identical, the curves of Pcsig) of 
the forward and the backward quench almost collapse onto the same curve (see Fig. lc), and the fidelity at t — 2(<?o — 9t)t~q is 
close to unity. 

We also plot the fidelity as a function of the quench time scale tq (see Fig. 2a). It can be seen that in both the impulse regime 
(tq <~ 10~ 3 ) and the adiabatic regime (tq >~ 10 2 ' 5 ), the fidelity is equal to unity. This agrees with our intuition. Meanwhile, 
in the intermediate regime, ~ 10~ 3 < tq <~ 10 2 5 , the fidelity oscillates rapidly (see Fig. 2a). When we plot the fidelity as 



a function of the quench time tq, instead of In tq, we found that there is a regular quasi-periodic oscillation (see Fig. 3). 
obtain an accurate expression of fidelity in the intermediate regime, 



We 
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it/2 
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k>0 



e -27TT Q sin 2 fegi^fc , 



27ttq sin k 
T 2 (l-iT Q sin 2 k) 



_g— WTQSin k e ~i(f> k 



(ID 



where 0fe = 2tq[(— gT + cosfc) 2 +sin k In a/4tq (— gr + cosfe) 2 ], andT(l— irgsin k) is the Gamma function (see Appendix 



A for details of the derivation). From Fig. 3 it can be seen that the analytical results agree with the numerical simulations, and 
that the fidelity oscillates quasi-periodically with the increase of tq as expected. This oscillation can actually also be observed 
in Fig. 1 (see Fig. lb' and Fig. lc'). 

Numerical simulations agree with the results obtained from the KZM very well, i.e., they account for three regimes that 
correspond to different tq. We are especially interested in the first regime - the adiabatic regime. From Fig. la and Fig. lc, it 
can be seen that in both the adiabatic regime and the impulse regime, the fidelity is close to unity. In the next subsection, we will 
introduce a method to eliminate the "degeneracy" of the adiabatic regime and the impulse regime. 




FIG. 3: Fidelity as a function of the time scale tq. Here the black dashed line represents the numerical results whereas the red solid line 
represents the analytical results (Eq. (Hit). It can be seen that the analytical results (Eq. i ll lb ) agree well with the numerical results except 
for the case tq — > 0, where the quench process enters the impulse regime. The length of the spin chain is TV = 50, and the delay time at the 
turnaround point is At = 0. 



C. Free evolution and decay of fidelity 



To distinguish the adiabatic and the impulse regime using quench echo, one can let the system evolve freely for some time 
at the turnaround point before quenching back. A study of the Landau-Zener problem with waiting at the minimum gap has 
been reported in Ref. 112511 . It was observed that the waiting influences the excitation probability. Similarly in our study the free 
evolution at the turnaround point leads to a decay in the fidelity in the impulse regime, but makes no difference in the adiabatic 
regime (see Figs, la'-lc' and Figs. la"-lc"). The reason is straightforward. Let us first consider the adiabatic regime. Because 
the system is always in its instantaneous ground state, the effect of the free evolution is simply a global phase factor, which does 
not affect the fidelity (see Figs, la', la", and Figs. 2b-2d). In the impulse regime, the wavefunction before the free evolution 
is the ground state of the initial Hamiltonian H(go = 10), and alternatively, a superposition of the excited and the ground 
states of H(gx = 0). The excited and ground states acquire different phase factors during the free evolution. Thus the wave 
function acquires relative phase factors in its components and is no longer the ground state of H{cjt = 0), but a superposition 
of its ground and excited states. Hence, in the impulse regime when one quenches the system back to the initial Hamiltonian 
H(g = 10), the system will no longer be in its ground state, but in a superposition of the ground state and the excited states. 
As a result, the fidelity is less than unity (see Fig. lc', Fig. lc", and Figs. 2b-2d). The length of time of the free evolution At 
also influences the fidelity. One can analytically calculate the fidelity as a function of the time of free evolution At: 



\-(g = +ao)) k 



\+(9T)) k +cos 



■(9t)) h , 



where 9k = arctan( 



cos k—gT • 



After free evolution for At, the wave function becomes 



sin 



,-iA kM At \ + {gT))k + 



COS 



JA k (g T )At 



<9T)) k . 



The fidelity can then be calculated as 



-r-r / sin 2 /csin 2 [A A: ( 5T )At] 

fc>0 V 



2<?t cos k + cjj, 



(12) 



Note that for a fixed chain size, the value of At needed to scramble all the relevant phases is relevant to the range of the spectrum 
of the system or the size of gap A^gx) of different k at the turnaround point gx- When A^g?) is very small, i.e., the energy 
spectrum of the system is concentrated within a very small energy range, one needs to wait for a long time in order to scramble 
all the relevant phases: At is inversely proportional to the energy scale J of A& (<7t). For a spin chain of N = 50, when the time 
of free evolution is very short, e.g., At = 0.1, there is a pronounced decay in the fidelity in the impulse regime (see Fig. 2b). The 
analytical result gives F w 0.882, which agrees with the numerical result. The fidelity decreases with the increase of time of the 
free evolution. The fidelity decays to 0.002 when At = 0.7 (see Fig. 2d). Hence the quench echo with a free evolution at the 
turnaround point can distinguish the adiabatic and the impulse regime. Our numerical results confirm our theoretical predictions. 

IV. BEYOND THE LINEAR QUENCH 

In the above discussion, we focused on the linear quench. One may repeat the above process with different tq until one finds 
the smallest Tq, under which the process is sufficiently adiabatic, for example F > 0.9. Nevertheless, the linear quench with 
Tq obtained above may waste a lot of time. The reason is obvious: in different regions of the parameter g, the energy gaps are 
different. According to KZM, different energy gaps correspond to different relaxation time r. For a linear quench protocol, we 
are treating the whole range of the parameter uniformly, and the relaxation time is determined by the global minimal energy gap. 
Thus, we waste a lot of time. Usually we want to ensure that the process not only nearly adiabatic but also as fast as possible. In 
the following we will consider nonlinear quench. 

A. Adjusting quench rate to the instantaneous gap 

An improved scheme is to divide the whole range of the parameter into many, e.g., M, parts with equal length (go — gr)/M, 
and then apply the above linear quench protocol to these ranges separately to find the uniformly adiabatic quench for each range 
Tq, i = 1, 2, • ■ ■ , M. We can also use the KZM to find a uniformly adiabatic quench protocol. From the discussion in Section 
III. A we know that the transition time scale is given by the absolute value ofA(g(t))/4:A(g(t)). Meanwhile, the relaxation time 
scale is given by 1/A(g(t)). When the former is many times larger than the latter, the process should be uniformly adiabatic. 
That is, when the parameter g(t) satisfies the relation 



A( 5 (t)) 



£A(ff(t)) 



A(flf(t)) 



(13) 



where 7 is a constant many times larger than unity, e.g., 7 = 10, the process is uniformly adiabatic in the sense that the ratio of 
two time scales remains a constant. Such a quench scheme is better than the linear quench. The solution to the above ordinary 
differential equation is 



A(g(t)) 



1 



T-t + c' 



(14) 



where qp corresponds to the sign of A/A on the left-hand-side of Eq. ( TT3l ) being positive or negative, and c is a constant of 
integration. For simplicity c can be chosen such that at t = A in Eq. ( TT4b is the minimal gap. Now, we know exactly the 
energy gap as a function of the controlling parameter (see Fig. 4a) J26I1 



A( 5 (t)) =2J A /1 - 2g(t) cos(-) + g*(t). 



(15) 



Therefore, c can be determined by g(t = 0) = cos (ir/N). Combining Eqs. (fl4l > and (Tl~5T >. we find the following uniformly 
adiabatic quench protocol (see Fig. 4b) 
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2 J sin (?r /AT) 
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(7) 2 



4 J 2 -t 



r.(- 



2 J sin O /N) 



7 



2Jsin(7r/AT) 



(0 < t < 



<t<0) 



2Jsin(7r/iV) 



(16) 



It can be seen that the time required for the whole process (quenching the controlling parameter from g 
by 



Ai K z = 



7 77V 

Jsin(7r/iV) ~ Jn ' 



to g = oo ) is given 



(17) 



or Aixz = a 7 » which is proportional to the chain size N and the ratio 7, and inversely proportional to the minimum energy 



gap A min = 2 jsin -^ ss 2Jj,, In linear quench the minimal time required for the adiabatic evolution grows with the system 
size like TV 2 12411 . The quench scheme of Eq. ( TToT ) is obviously better. This agrees with previous studies that "non-linear" quench 
can improve the adiabaticity (minimize excitation) IU8l420ll . The energy gap and the protocols for uniformly adiabatic quench 
(Eq. ( fToT l) are shown in Fig. 4. 





g(t) 



FIG. 4: (a) Energy gap as a function of the controlling parameter g for a finite size chain. Here the solid line indicates the spin chain of 
N = 50 and the dashed line indicates the gap of an infinite chain, (b) The uniformly adiabatic quench protocol associated with the KZM 
criterion (solid) dl6b and associated with RC criterion (dashed) d!9t . Here the ratio 7 = 2, and the energy scale J = 1/2. We have chosen the 
condition g(t — 0) = cos (n/N). 



Note that the criterion for uniformly adiabatic evolution derived from the KZM (Eq. ( fT3] )) is similar, but not identical, to the 
criterion proposed by Roland and Cerf that was derived from the quantum adiabatic theorem (See Eq. (17) of reference 12711 ). 
In the RC model the energy gap is inversely proportional to \/N, and the minimum time required is proportional to y/~N. But 
in the Ising chain, the energy gap is inversely proportional to N, and the minimum time required is proportional to N. It can be 
proved that if one uses Roland and Cerf's criterion to evaluate the minimum time required for the uniformly adiabatic evolution, 
the minimum time is also proportional to N. 

It is interesting to compare the two criteria for uniformly adiabatic evolution in the Ising chain. In the following we will 
first solve the equation of the quench protocol for uniformly adiabatic evolution .9rc(£) associated with the Roland and Cerf's 
criterion and then simulate the dynamic evolution of the Ising chain with both gnc(t) an d 9Kz(t)- We will fix the time of quench 
process, and compare the fidelity of the two protocols. The Roland and Cerf's criterion (see Eq. (17) of Ref. 112711 ) is 



dt 



\g(t) - g c 



1 



A 2 ( 5 W), 



(18) 



where 7' is the ratio between the two time scale. Obviously, when N —$■ 00, A(g) = \g(t) — g c \ is valid for arbitrary g. In this 
respect, the two criteria, Eq. (TT~3b and Eq. ( fT8l , are equivalent. Nevertheless, when N is finite, the two criteria differ slightly 
because the gap A(g) deviates from \g(t) — g c \ near the critical point (see Eq. ( TT3T > and Fig. 4a). As a result there is a small 
discrepancy in the quench protocols <?rc(£) and <?kz(£) associated with two criteria, especially when g(t) is close to g c . 
By substituting Eq. (fT~5b into Eq. ( fT8l . we obtain the quench protocol: 



9Rc{t) = cos( — ) +sin(— )tan 



2Jsin(f) i 






4.7 



4.7 



(19) 



We plot the solution <7rc(X) along with gKz(i) in Fig- 4b. There is a "kink" at the anti-cross point of the energy levels in gKz(t) 
associated with the KZM criterion, but there is none in gnc(t) associated with the RC criterion (see Fig. 4b and the inset). 
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Although there is a singularity (divergent time derivative of g) in g^z (t) at t = 0, the time interval of this region is vanishingly 
small. As a result the total change in g in this singular region is very small (See Eq. 4b), and the eigenstates of H(t) do not 
change significantly within it. This is in the same spirit as quantum fidelity {3Qi], where fidelity susceptibility diverges at quantum 
critical point, but the fidelity is nonzero indicating that the ground state does not changes significantly. Hence, the 'kink' at t = 
will not lead to a lot of excitations. Our simulation verifies this point. Similar to Eq. 1171 . we obtain the time required for the 
uniformly adiabatic evolution 



At 



RC 



7'JV 
~2J"' 



(20) 



Comparing Eq. ( TT7T ) and Eq. d20l i. we find that when 7' = -7, the time required for two criteria are equal. In the following we 
will simulate the dynamics of the Ising chain under the two quench protocols: Eq. ( fT6] l and Eq. ( fT9l ). Substituting Eq. ( fT6] l and 
Eq. ( fT9l into Eq. ©, one obtains the instantaneous fidelity as a function of the time F = \j3(t)\ associated with two criteria. 
We plot the fidelity as a function of the time in Fig. 5. When one chooses a different initial condition, the fidelity as a function of 
the time differs a lot. In the left panel of Fig. 5, we plot the fidelity as a function of the time quenching from t = ~ 21 sinful N) 
(g = —00). The fidelity associated with the Roland and Cerf criterion decays when the system is near the anti-crossing point, 
and then revives. But the fidelity associated with the KZM does not change much, and remains close to unity all the time. At 
t = 2jsin(Ti/N) (9 = °°)' tne fidelity associated with the Roland and Cerf 's criterion is a bit higher than that associated with the 



KZM. In the right panel of Fig. 5, we plot the fidelity as a function of the time starting from t = — -k 2/ S infa/AO • ^ n ^ s case 
the fidelity associated with the Roland and Cerf's criterion oscillates rapidly and finally reaches a stable value around 0.85. By 
contrast the fidelity associated with the the KZM does not oscillate and remains very close to unity. From the above facts, we 
conclude that in some cases, the Roland and Cerf's criterion is better than the KZM criterion, but in some other cases, it is worse. 
Hence, we cannot say which criterion is definitely better, but the KZM provides new insights into the conditions for uniformly 
adiabatic evolution. 

One is tempted at this point to undertake a variational study in search of optimal quenches. While such a study is beyond the 
scope of this paper, we note that in practical applications (e.g., adiabatic quantum computing) optimization would involve not 
just varying rate, but (as it was done in Fig. 5) also the starting and final points of the quenches can be brought closer to the 
"critical point". Resulting errors can be detected and the correct result can be ascertained by repeating the computation many 
times. 



Fidelity 
-14). 



Fidelity 




-300 -200 -100 100 200 300 -100 



FIG. 5: Comparison of the Roland and Cerf criterion and the KZM criterion. The horizontal axis indicates the time, and the vertical axis 
depicts the instantaneous fidelity. The solid (dashed) line describes the fidelity as a function of the time t associated with the Roland and Cerf 
(KZM) criterion. Left panel: The evolution starts from the ground state of t — —jj^ (g = —00). Right panel: The evolution starts from the 
ground state of t 



1jN_ 
' 3 2Jtv ■ 



B. Gauging the distance from the adiabatic quench 



Fig. 2 indicates that when the time scale of the quench tq is in the range tq <G (10 _1 , 1), the fidelity is almost equal to zero. 
However, this does not reveal how far the quench is from the adiabaticity. For example when one out of many (N = 50 in our 
numerical simulation) modes get excited, the fidelity will decay to nearly zero due to the orthogonality of one mode. But, in a 
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sense, the system is still close to the ground state, as all but one excited state are empty. In this sense, the fidelity is not a good 
criterion for measuring how far away the quench is from the adiabatic evolution. 

A better gauge of the distance of a quench process from the adiabaticity may be obtained using other variables, such as the 
magnetization per site along the direction of the external magnetic field 



N N 



(21) 



fc=i 



When the magnetic field is large, the ground state corresponds torn = 1. We plot the final magnetization as a function of tq 
in Fig. 6. In the range of tq <G (10 _1 , 1), the fidelity is vanishingly small, but the magnetization per site is still large. This 
indicates that the system is not very far away from the instantaneous ground state. Moreover, when one delays for some time 
at the turnaround point, the magnetization per site of the impulse regime will decrease, but that of the adiabatic regime will not 
(see Fig. 6). This is similar to the fidelity and agrees with our intuition. Last but not least, the magnetization is experimentally 
easier to the measure than the fidelity, and it has been used as a tool to study the adiabaticity of quantum dynamics in Ref. 12911 . 
One can also use the kink density | 2,(1 — (cjCTj+i)) [240 as a measure of the distance of the system from the adiabaticity. 
The relation between the fidelity and the density of defects has been studied in Ref. 113111 . Other variables, such as the residual 
energy 0211 . can be also used to gauge the distance from the adiabatic quench. Such obvious measures of how far the quench is 
from the adiabaticity work well in the one-dimensional Ising model, but finding their useful analogues in other situations (e.g., 
adiabatic quantum computing) may not be easy. 
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FIG. 6: Final magnetization of the Ising chain as a function of the time scale of the quench tq. All the parameters are the same as those 
in Fig. 2(a). One can see that the spin chain is not very far away from equilibrium except when the time scale of the quench is in the range 
tq e (1CT 1 , 1). Left panel: delay time At = 0; Right panel: delay time At = 0.7. 



V. SUMMARY AND CONCLUSION 

We have proposed a strategy to test the adiabaticity without knowing either the eigenstates or eigenenergies of the Hamiltonian. 
Instead of having to find the gap of the Hamiltonian, and then using the quantum adiabatic approximation to evaluate the 
adiabaticity, one can use a quench echo to evaluate the adiabaticity of an evolution. The underlying mechanism is that when the 
time scale of the quench is large in comparison with the inverse of the energy gap, both the forward and the backward evolutions 
are adiabatic. As a result, the fidelity of the initial state and the final state is close to unity. Otherwise, the evolution is not 
adiabatic, and the fidelity is less than unity. The method for testing the adiabaticity of an evolution presented in this paper is 
universally valid. It does not depend the model or the validity of conditions for adiabatic approximation. We further proposed 
a method for finding the uniformly adiabatic quench protocol based on the KZM, and discussed the problem of gauging how 
non-adiabatic is a quench. Given the importance of the adiabaticity in various applications, we believe that our results will be 
broadly applicable, and may be useful in experimental applications. 
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Appendix A: fidelity in the intermediate regime 

From Refs. 11241,12811 . we know that in the wave vector k representation, the Schrodinger equation for the forward quench can 
be rewritten as Landau-Zener type equations (see Eq. (|6]l for a comparison): 



. d 

l dF 



v k (f) 

Ufc(t') 



4- i 
1 -4- 



Vk(t') 

Uk(t') 



(Al) 



where t' = 4tq smk[— g(t) + cos k] and t'q = 4tq sin 2 k. This equation can be solved in terms of Weber functions. The initial 
conditions are Vk(t' = — oo) = and Uk(t' = — oo) = 1. The solution for this equation is 112811 



uie(t) «e~ 7rTQSin2fc e ?TQ "~ 9+cosA; ' >2+sin2 k ln V iT Q (-9+cos k) 2 ] 

\2TTTnsm 2 k x T05i „2 fc r , s9 . , _ ( A2 ) 

Vk U\ ~ V W e ^ e - iT Q l(-9+c°s fc)-+sin 2 k ln ^4tq (- 9 +cos fc)»] _ 

T(l — itq sin 2 A;) 

Similarly, we can obtain the solution for the quench echo. Combining the forward and the quench echo process, we find the 
solution of the fidelity (fTTT i. It is worth pointing out that the above solution is only good for the turnaround point far away from 
the critical point gx <C 1. For example gx = 0.5. 
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